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Abstract 

Simple analytically solvable models are proposed exhibiting l/f spectrum in wide 
range of frequency. The signals of the models consist of pulses (point process) which 
interevent times fluctuate about some average value, obeying an autoregressive pro- 
cess with very small damping. The power spectrum of the process can be expressed 
by the Hooge formula. The proposed models reveal possible origin of 1// noise, i.e., 
random increments of the time intervals between pulses or interevent time of the 
process (Brownian motion in the time axis). 

1 Introduction 

Widespread occurrence of signals exhibiting power spectral density with l/f behavior 
suggests that a general mathematical explanation of such an effect might exist. However, 
physical models of l/f noise in some physical systems are usually very specialized and 
they do not explain the omnipresence of the processes with 1/ f s spectrum. Mathematical 
algorithms and models for the generation of the processes with l/f noise can not, as a 
rule, be solved analytically and they do not reveal the origin as well as the necessary and 
sufficient conditions for the appearance of l/f type fluctuations. 

History of the progress in different areas of physics indicates to the crucial influence 
of simple models on the advancement of the understanding of the main points of the new 
phenomena. We note here only the decisive influence of the Bohr model of hydrogen atom 
on the development of the quantum theory, the role of the Lorenz model as well as the 
logistic and standard (Chirikov) maps for understanding of the deterministic chaos and 
the quantum kicked rotator for the revealing the quantum localization of classical chaos. 

Here we present simple models of l/f noise which may essentially influence on the 
understanding of the origin and main properties of the effect. Our models are the result 
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of a search for necessary and sufficient conditions for the appearance of 1/f fluctuations 
in simple systems affected by the random external perturbations, which where initiated 
in jH [2] and originated from the observation of the transition from chaotic to nonchaotic 
behavior in the ensemble of randomly driven systems [3J E] ■ 



2 Simple model 



The simplest version of our model corresponds to one particle moving along some orbit. 
The period of this motion fluctuates (due to external random perturbations of the system's 
parameters) about some average value f . So, a sequence of the transit times when 
the particle crosses some point of the orbit is described by the iterative equations 



tk-l + T ki 

Tfc-i -7(^-1 -f) +ae k . 



(1) 



Here 7 is the coefficient of the relaxation of the period to the average value f , {e fc } denotes 
the sequence of random variables with zero expectation and unit variance and a is the 
standard deviation of the noise. Due to the contribution of a large number of random 
variables to the transit times, model (1) represents a long-memory random process. 

Power spectral density S (f) of the signal or current of the model (1), I (t) = a J2k $ (t ~ tk) 
(with a being a contribution to the signal of one pulse or contribution to the current of 
one particle when it crosses the section of observation), may be calculated according to 
equation 




S(f)= lim — 



E 



-i2nft k 



(2) 



where T is the whole observation time interval, k m i n and /c max are minimal and maximal 
values of index k in the interval of observation and the brackets (...) denote the averaging 
over the realizations of the process. Eq. (2) may be rewritten in the form 



S(f) 




i2irfA(k;q) 



(3) 



where A (k; q) = tk +q — tj, is the difference of transit times tk +q and t/.. 

At k ^> 7 _1 we have the stationary process the expectation (rfc) = r and the 
variance = (rf) — (r^) 2 = a 2 /2'y of the recurrence time do not depend on the time 
while A (k; q) is a normally distributed random variable with the expectation /j,a {q) = 
(A (00; q)) = qf and the variance a\ (q) = (A (00; q)j — (A (00; q)) 2 expressed as 



*i(g) 



(l-(l-7f) 
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For the normal distribution of A (k; q) = A (q) Eq. (3) yields 

S(f) =2 /a^e i2 ^ (9) -^ 2/2<J i ((3) 



(4) 



(5) 
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where I = lim a (fc max — k min + I) /T = a/r is the average current. 

T — >oo 

Substitution of the expansions of the variance (4) at \q\ <C 7 _1 in powers of 7 \q\, 

°\ (<?) = ^r? 2 , 
into Eq. (5) yields 1 / /-like power spectrum, 

s(f) = i 2 Y> ( 6 ) 

for sufficiently small parameters a and 7 in any desirably wide range of frequencies, f\ = 
7/71-0Y < / < f 2 — l/na T . Here an is a dimensionless constant (the Hooge parameter), 

a H = ^=Ke- K \ K=-J=-. (7) 
V 71 " \j2o T 



We see that the power of 1/ f noise except the squared average current I 2 depends strongly 
on the ratio of the average recurrence time f to the standard deviation of the recurrence 
time a T . 

Therefore, the process (1) containing only one relaxation time 7" 1 can for sufficiently 
small damping 7 produce an exact 1 / /-like spectrum in wide range of frequency {fx, f 2)1 
with /2//1 ~ 7 _1 . 

The model is free from the unphysical divergency of the spectrum at / — > 0. So, using 
for / < /1 an expansion of expression (4) at \q\ ^> 7^, o\(q) = (cr/7) 2 (|g| — I/7), we 
obtain from Eq. (5) the Lorentzian power spectrum [Hj 



S(f) = I 



2 4T re ; 

1 + tLu 2 ' 



'rel L 

Here u = 2nf and T re i — D t — a 2 /2f , y 2 is the "diffusion" coefficient of the time For 
/ ^ fo — f^/ncr 2 = l/2nr re i we have the white noise 

S(f) = P {2a 2 1 f^) , 
/ < min {fx, f = T-f 2 /na 2 } . 

This is in agreement with the statement |& that the power spectrum of any pulse sequence 
is white at low enough frequencies. 

Equations (6)-(8) describe quite well the power spectrum of the random process (1) for 
perturbation by the Gaussian white noise and even for perturbations by the non-Gaussian 
sequence of random impacts {e^} in Eq. (1) (see illustrative examples in Refs. [S] and jjj). 



3 Generalizations and numerical analysis 

This simple exactly solvable model can easily be generalized in different directions: for 
large number of particles moving in similar orbits with coherent (identical for all particles) 
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or independent (uncorrelated for different particles) fluctuations of the periods, for non- 
Gaussian or continuous perturbations of the systems' parameters and for spatially extended 
systems. So, when an ensemble of N c particles moves on closed orbits and the period of 
each particle fluctuates independently (due to the perturbations by uncorrelated sequences 
of random variables {e k }, different for each particle v) the power spectral density of the 
collective current / of all particles can be calculated by the above method [5 too and it is 
expressed as the Hooge formula |8>, 9 a 

S{f) = I 2 ^- f N c = Vn c . (6a) 

This expression with the factor 1/N C is in agreement with the Hooge statements 
that summation of spectra is only allowed if the processes contributing to the spectrum 
are isolated from each other and that only isolated traps yield 1/ / spectrum. 

It should be noted that in many cases the intensity of signals or currents can be ex- 
pressed as a sequence of the pulse occurrence times {tk}, i.e., as I (t) = aJ2k$ (t ~ tk)- 
This expression represents exactly the flow of identical point objects: cars, electrons, pho- 
tons and so on. More generally, instead of the Dirac delta function one should introduce 
time dependent pulse amplitudes A k (t — tk). The low frequency power spectral density, 
however, depends weakly on the shapes of the pulses P, while fluctuations of the pulses 
amplitudes result, as a rule, in white or Lorentzian, but not 1//, noise. Model (1) in such 
cases represents fluctuations of the averaged interevent time Tk between the subsequent 
occurrence times of the pulses. 

The model may also be generalized for the nonlinear relaxation of the interevent time 
Tk- In such Eq. (1) can be written in the form 

= tk-i + n, 

rfvfa-i) , (9) 

Here the function V(tjj.) represents the effective "potential well" for the Brownian motion 
of the interevent time 7*. The steady state distribution density of the interevent time Tk 
generated by Eq. (9) is of the form 



V?(7*) = Cexp 



2V(r k ) 
a 2 



(10) 



where a constant C may be obtained from the normalization. For the power-law "potential 
well" 

v(n) = ^(T k -ff n (li) 

with integer n we have a generalization of Eqs. (1) 

tk — tk-l + Tk, n „x 

/ -\2n-l , K 1 -*) 

Tk = T k -1 - in (T k -1 -T) + (T6 k . 
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The steady state distribution density of the interevent time r k in such a case is 

%W- 2 r(l + l/2n) ' 1 ' 

For sufficiently large n ^> 1 Eqs. (12) and (13) represent Brownian motion of the interevent 
time r k in almost rectangle "potential well" (square-well potential) restricting movement 
of r k mostly in the interval (f — h, f + h) with h ~ (a 2 /7) 1//2n . At h < f such a restriction 
prevents from the emergence of the negative interevent times r k and, consequently, from 
the clustering of the particles or of the signal pulses. 

We can evaluate the power spectrum of the processes (9) and (12) as well [SI CI]- The 
power spectral density (3) may be rewritten in the form 

S (f) = 21 a {J2 e i2nfTk{q)q ^ (14) 

where the occurrence times t k+q and t k difference A (k; q) is expressed as 

k+q 

A (k; q) = t k+q -t k = T k = r k (q)q,q>0 

l=k+l 

with r k (q) = (tk+q — tk)/q being the average interevent time in the time interval (t k ,t k+q ) 
and the brackets in Eq. (14) denote averaging over time (index k) and over realizations of 
the process. At 2^fr k (q) <C 1 the summation in (14) may be replaced by the integration 



-t-cxj 

S {f) = 21 a J (e i2nfTkq ) dq. 



Here the averaging over k and over realizations of the process coincides with the averaging 
over the distribution of the interevent times r k , i.e., 



-oo 



e i2nfq^ = J e i2nfqr^ ^ ^ = ^ ^ f q) 
-co 

where ip T (r k ) is the distribution density of the interevent times r k and \ T (i?) is the char- 
acteristic function of this distribution. 

Taking into account the property of the characteristic function 

Xr (#) d$ = 2vn/v (0) 

we obtain the final expression for the power spectrum 

S (f) = 2Pf (0) //. (15) 
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For Gaussian distribution of the interevent times T k 

ij T (0) = exp (-f 2 /2al) j^a T 

and according to Eq. (15) we recover expressions (6) and (7) obtained from the analysis 
of the dynamical process. 

Therefore, according to Eqs. (6), (7) and (15) the dimensionless parameter an = 
2fip T (0) is proportional to the distribution density ip T (0) of the interevent time Tk in 
the point r k = 0, i.e., to the probability of the clustering of the signal pulses. The pulse 
clustering results in the large variance of the signal - the condition necessary for appearance 
of stationary 1/f noise in the wide range of frequency. 

It should be noted, however, that Eq. (15) represent an idealized 1/f noise. The 
real systems have finite relaxation time and, therefore, expression of the noise intensity in 
the form (14)-(15) is valid only for / > (27rr re ;) _1 with r re / being the relaxation time of 
the interevent time's fluctuations. On the other hand, due to the deviations from the 
approximation t k+q — tk = T^q at large q, for sufficiently low frequency we can obtain the 
finite intensity of 1/f 5 (5 ~ 1) noise even in the case ip (0) = but for the signals with 
fluctuations resulting in the dense concentrations of the pulse occurrence times tk- 

We can generate, of course, the stationary time series of the occurrence times tk also for 
other restrictions for the interevent time Tfc, e.g., with the reflecting boundary conditions 
at some values r k = r min and r k = r max . The process like (1) with the reflecting condition 
for T k at r min = r s may also be generated by the recurrent equations 

tk = tk-l + r s + Tk ,-.„>. 
Tk = \r k -i - l(rk-i - f) + ae k \ - 

Numerical analysis of the models like (1), (9), (12) and (16) as well as with other 
restrictions for the interevent time Tk shows that power spectrum of the current is 1/ /-like 
in large interval of frequency only when the distribution density of the interevent times T k 
in the point Tk — is nonzero, i.e., ipri^k — 0) 7^ 0. For models with ip T (0) = or Vv(0) 
very close to zero we observe in numerical simulations the power spectrum S (/) oc l// 3//2 
(see Ref. jT] and Fig. 1). 

4 Conclusions 

Simple analytically solvable models of 1/f noise are proposed. The models reveal main 
features, parameter dependences and possible origin of 1/f noise, i.e., random increments 
of the time intervals between the pulses or interevent times of the elementary events of the 
process. The conclusion that 1/f 6 noise with S ~ 1 may result from the clustering of the 
signal pulses, particles or elementary events can be drawn from the analysis of the simple, 
exactly solvable models. The mechanism of the clustering depends on the system. 
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Fig. 1. Power spectrum vs frequency of the current generated by Eqs. (l)-(3) and (16) with 
parameters f = 100, a = 1 and with the Gaussian distribution of the random increments 
The sinuous curves represent the results of numerical simulations: (a) and (b) according to Eq. 
(1) with 7 = and with reflecting boundary conditions for 77 at f — h and at r + h for different 
values of h; (c) and (d) according to Eq. (16) with 7 = 0.0001 and different values of t s . 
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